Robust entanglement of a micromechanical resonator with output optical fields 



OO 
O 

o 

(N 
C 

(N 



^ ■ 
-t— > 



> 
in 

o 

(N 

o 

oo 
o 



X 



C. Genes, A. Mari, P. Tombesi, and D. Vitali 

Dipartimento di Fisica, Universita di Camerino, 1-62032 Camerino (MC), Italy 

We perform an analysis of the optomechanical entanglement between the experimentally de- 
tectable output field of an optical cavity and a vibrating cavity end-mirror. We show that by a 
proper choice of the readout (mainly by a proper choice of detection bandwidth) one can not only 
detect the already predicted intracavity entanglement but also optimize and increase it. This en- 
tanglement is explained as being generated by a scattering process owing to which strong quantum 
correlations between the mirror and the optical Stokes sideband are created. All-optical entangle- 
ment between scattered sidebands is also predicted and it is shown that the mechanical resonator 
and the two sideband modes form a fully tripartite-entangled system capable of providing practicable 
and robust solutions for continuous variable quantum communication protocols. 

PACS numbers: 03.67.Mn, 85.85.+j,42.50.Wk,42.50.Lc 



I. INTRODUCTION 

Mechanical resonators at the micro- and nano-meter 
scale are now widely employed in the high-sensitive de- 
tection of mass and forces [3, S Q ■ The recent improve- 
ments in the nanofabrication techniques suggest that in 
the near future these devices will reach the regime in 
which their sensitivity will be limited by the ultimate 
quantum limits set by the Heisenberg principle, as first 
suggested in the context of the detection of gravitational 
waves by the pioneering work of Braginsky and coworkers 

a. 

The experimental demonstration of genuine quantum 
states of macroscopic mechanical resonators with a mass 
in the nanogram-milligram range will represent an im- 
portant step not only for the high-sensitive detection 
of displacements and forces, but also for the founda- 
tions of physics. It would represent, in fact, a remark- 
able signature of the quantum behavior of a macro- 
scopic object, allowing to shed further light onto the 
quantum-classical boundar y [j l . Si gni ficant experimental 
[I0,B[I[lMM^ and 
theoreticallllllinilWtla 113 HE M SWefForts 
are currently devoted to cooling such microresonators to 
their quantum ground state. 

However, the generation of other examples of quan- 
tum states of a micro-mechanical resonator has been 
also considered recently. The most relevant examples 
are given by squeezed and entangled states. Squeezed 
states of nano-mechanical resonators [HJ are potentially 
useful for surpassing the standard quantum limit for po- 
sition and force detection 0], and could be generated 
in different ways, using either the coupling with a qubit 
[33| , or measurement and feedback schemes [25|, [34[ . En- 
tanglement is instead the characteristic element of quan- 
tum theory, because it is responsible for correlations be- 
tween observables that cannot be understood on the ba- 
sis of local realistic theories 35] . For this reason, there 
has been an increasing interest in establishing the con- 
ditions under which entanglement between macroscopic 
objects can arise. Relevant experimental demonstration 



in this directions are given by the entanglement between 
collective spins of atomic ensembles [36| . and between 
Josephson-junction qubits [37|]- Then, starting from the 
proposal of Ref. [38| in which two mirrors of a ring cav- 
ity are entangled by the radiation pressure of the cav- 
ity mode, many proposals involved nano- and micro- 
mechanical resonators, eventually entangled with other 
systems. One could entangle a nanomechanical oscillator 
with a Cooper-pair box [3{|, while Ref. [4(| studied how 
to entangle an array of nanomechanical oscillators. Fur- 
ther proposals suggested to entangle two charge qubits 
[4l| or two Josephson junctions [42| via nanomechani- 
cal resonators, or to entangle two nanomechanical res- 
onators via trapped ions (43[, Cooper pair boxes [3], 
or dc-SQUIDS More recently, schemes for entan- 

gling a superconducting coplanar waveguide field with a 
nanomechanical resonator, either via a Cooper pair box 
within the waveguide [46| . or via direct capacitive cou- 
pling [47|, have been proposed. 

After Ref. [11], other optomechanical systems have 
been proposed for entangling optical and/or mechani- 
cal modes by means of the radiation pressure interac- 
tion. Ref. 48] considered two mirrors of two different 
cavities illuminated with entangled light beams, while 
Refs. [i^, [5(| HH H2] considered different examples of 
double-cavity systems in which entanglement either be- 
tween different mechanical modes, or between a cavity 
mode and a vibrational mode of a cavity mirror have been 
studied. Refs. [H, [54| considered the simplest scheme 
capable of generating stationary optomechanical entan- 
glement, i.e., a single Fabry-Perot cavity either with one 
[53| . or both (54[, movable mirrors. 

Here we shall reconsider the Fabry-Perot model of 
Ref. [53|], which is remarkable for its simplicity and ro- 
bustness against temperature of the resulting entangle- 
ment, and extend its study in various directions. In fact, 
entangled optomechanical systems could be profitably 
used for the realization of quantum communication net- 
works, in which the mechanical modes play the role of lo- 
cal nodes where quantum information can be stored and 
retrieved, and optical modes carry this information be- 
tween the nodes. Refs. [H, [H, H3] proposed a scheme of 
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this kind, based on free-space light modes scattered by 
a single reflecting mirror, which could allow the imple- 
mentation of continuous variable (CV) quantum telepor- 
tation [55( | , q uantum telecloning [56[, and entanglement 
swapping [|57[. Therefore, any quantum communication 
application involves traveling output modes rather than 
intracavity ones, and it is important to study how the op- 
tomechanical entanglement generated within the cavity 
is transferred to the output field. Furthermore, by con- 
sidering the output field, one can adopt a multiplexing 
approach because, by means of spectral filters, one can 
always select many different traveling output modes orig- 
inating from a single intracavity mode (see Fig. 1). One 
can therefore manipulate a multipartite system, eventu- 
ally possessing multipartite entanglement. We shall de- 
velop a general theory showing how the entanglement be- 
tween the mechanical resonator and optical output modes 
can be properly defined and calculated. 

We shall see that, together with its output field, 
the single Fabry-Perot cavity system of Ref. [53[ rep- 
resents the "cavity version" of the free-space scheme of 
Refs. [E^,|5t|. In fact, as it happens in this latter scheme, 
all the relevant dynamics induced by radiation pressure 
interaction is carried by the two output modes corre- 
sponding to the first Stokes and anti-Stokes sidebands of 
the driving laser. In particular, the optomechanical en- 
tanglement with the intracavity mode is optimally trans- 
ferred to the output Stokes sideband mode, which is how- 
ever robustly entangled also with the anti-Stokes output 
mode. We shall see that the present Fabry-Perot cavity 
system is preferable with respect to the free space model 
of Refs. [H, Hll, because entanglement is achievable in a 
much more accessible experimental parameter region. 

The outline of the paper is as follows. Sec. II gives 
a general description of the dynamics by means of the 
Quantum Langevin Equations (QLE), Sec. Ill analyzes 
in detail the entanglement between the mechanical mode 
and the intracavity mode, while in Sec. IV we describe a 
general theory on how a number of independent optical 
modes can be selected and defined, and their entangle- 
ment properties calculated. Sec. V is for concluding 
remarks. 



II. SYSTEM DYNAMICS 

We consider a driven optical cavity coupled by radia- 
tion pressure to a micromechanical oscillator. The typ- 
ical experimental configuration is a Fabry-Perot cavity 
with one mirror much lighter than the other (see e.g. 
[1, [111 [Sill, [IS [13]), but our treatment applies to other 
configurations such as the silica toroidal microcavity of 
Refs. [14|, |l9|, [58| . Radiation pressure couples each cavity 
mode with many vibrational normal modes of the mov- 
able mirror. However, by choosing the detection band- 
width so that only an isolated mechanical resonance sig- 
nificantly contributes to the detected signal, one can re- 
strict to a single mechanical oscillator, since inter-mode 
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FIG. 1: Scheme of the cavity, which is driven by a laser and 
has a vibrating mirror. With appropriate filters one can select 
N independent modes from the cavity output field. 



coupling due to mechanical nonlinearities are typically 
negligible (see also [59| for a more general treatment). 
The Hamiltonian of the system reads [(30[ 

H = tvjj c a) a + -hu> m (p 2 + q 2 ) — HGqo) aq (1) 



-ihE(a'e 



ae 



iu> t 



(2) 



The first term describes the energy of the cavity mode, 
with lowering operator a ([a, a'] — 1), cavity frequency 
lu c and decay rate n. The second term gives the en- 
ergy of the mechanical mode, modeled as harmonic os- 
cillator at frequency uj m and described by dimensionless 
position and momentum operators q and p {[q,p] = i)- 
The third term is the radiation-pressure coupling of rate 
Go = (uj c / L)y/h/ muj m , where m is the effective mass 
of the mechanical mode [6lJ , and L is an effective length 
that depends upon the cavity geometry: it coincides with 
the cavity length in the Fabry-Perot case, and with the 
toroid radius in the case of Refs. OH!]- The last term 
describes the input driving by a laser with frequency 
ujq, where E is related to the input laser power P by 
\E\ = y / 2PK/haj . One can adopt the single cavity mode 
description of Eq. @ as long as one drives only one cavity 
mode and the mechanical frequency cu m is much smaller 
than the cavity free spectral range FSR ~ c/L. In this 
case, scattering of photons from the driven mode into 
other cavity modes is negligible [6^ |. 

The dynamic is also determined by the fluctuation- 
dissipation processes affecting both the optical and the 
mechanical mode. They can be taken into account in a 
fully consistent way [6fj] by considering the following set 
of nonlinear QLE, written in the interaction picture with 
respect to favosa 



q = bJmP, 

p = -uj m q - j m P + Go<J a + £, 
a = —(K + iAo)a + iGoaq + E + 



(3a) 
(3b) 

■>i;a m , (3c) 



where Ao = lu c — ujq. The mechanical mode is affected by 
a viscous force with damping rate 7 m and by a Brownian 
stochastic force with zero mean value £, that obeys the 
correlation function [6(| [HI 



mm) 



7m / -iw(t-i') 

2tt 



coth 



2k B T J 



1 



(4) 
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where fcs is the Boltzmann constant and T is the temper- 
ature of the reservoir of the micromechanical oscillator. 
The Brownian noise £ (t) is a Gaussian quantum stochas- 
tic process and its non-Markovian nature (neither its cor- 
relation function nor its commutator are proportional to 
a Dirac delta) guarantees that the QLE of Eqs. ([3|) pre- 
serve the correct commutation relations between oper- 
ators during the time evolution [60]. The cavity mode 
amplitude instead decays at the rate k and is affected by 
the vacuum radiation input noise a m (t), whose correla- 
tion functions are given by 



and 



(a m (t)a m ^(t')) = [N(u D ) + 1] S(t - f)> ( 5 ) 



(a in ^(t)a in (t'))=N(uj c )6(t-t'), (6) 



where N(u> c ) = (cxp{huj c / ksT} — 1) 1 is the equilib- 
rium mean thermal photon number. At optical frequen- 
cies hwc/ksT 3> 1 and therefore N{u> c ) ~ 0, so that only 
the correlation function of Eq. ((SJ) is relevant. We shall 
neglect here technical noise sources, such as the ampli- 
tude and phase fluctuations of the driving laser. They 
can hinder the achievement of genuine quantum effects 
(see e.g. |19|). but they could be easily accounted for by 
introducing fluctuations of the modulus and of the phase 
of the driving parameter E of Eq. flS]) [H[ . 



A. Linearization around the classical steady state 
and stability analysis 

As shown in [53j| . significant optomechanical entan- 
glement is achieved when radiation pressure coupling is 
strong, which is realized when the intracavity field is very 
intense, i.e., for high-finesse cavities and enough driving 
power. In this limit (and if the system is stable) the 
system is characterized by a semiclassical steady state 
with the cavity mode in a coherent state with amplitude 
a s = E/(k + iA), and the micromechanical mirror dis- 
placed by q s = Go\a s \ 2 /uj m (see Refs. [13, HH HI] for 
details). The expression giving the intracavity ampli- 
tude a s is actually an implicit nonlinear equation for a s 
because 



A = An- 



(7) 



is the effective cavity detuning including the effect of the 
stationary radiation pressure. As shown in Refs. [30I.I53T]. 
when \a s \ 3> 1 the quantum dynamics of the fluctua- 
tions around the steady state is well described by lin- 
earizing the nonlinear QLE of Eqs. ([3|). Defining the 
cavity field fluctuation quadratures SX = (Sa + Sa^ ) j \[2 
and SY = (Sa — SaJ)/iv2, and the corresponding Hermi- 
tian input noise operators X m = (a m + a m, ^)/\/2 and 
yin - (a in - a in '^)/iy/2, the linearized QLE can be writ- 
ten in the following compact matrix form [53| 



where u T (t) = (Sq(t), Sp(t), SX(t), 5Y(t)) T ( T denotes 
the transposition) is the vector of CV fluctuation op- 
erators , n T (t) = {0,£(t),V^X m (t),V2^Y m (t)) T the 
corresponding vector of noises and A the matrix 



A = 


















G 











— K 
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G 





-A 


— K 



(9) 



where 



G = G a s V2 = 



2bJ r 



Pk 



L y mu m uo (k 2 + A 2 ) ' 



(10) 



is the effective optomechanical coupling (we have chosen 
the phase reference so that a s is real and positive) . When 
a s ^> 1, one has G ^ Go, and therefore the generation 
of significant optomechanical entanglement is facilitated 
in this linearized regime. 

The formal solution of Eq. © is u(t) = M(t)u(0) + 

J* dsM(s)n(t—s), where M(t) = exp{At}. The system is 
stable and reaches its steady state for t — > oo when all the 
eigenvalues of A have negative real parts so that M(oo) = 
0. The stability conditions can be derived by applying 
the Routh-Hurwitz criterion [67j|, yielding the following 
two nontrivial conditions on the system parameters, 



si 



2j m K | k 2 + (w m - A) 2 k 2 + (u m + A) 



+7m [(7m + 2k) (k 2 + A 2 ) + 2nuj 2 m ] } 
+Aiv m G 2 ( 7m + 2k) 2 > 0, 



S2 



(k 2 



G A > 0. 



(11a) 
(lib) 



which will be considered to be satisfied from now on. No- 
tice that when A > (laser red-detuned with respect to 
the cavity) the first condition is always satisfied and only 
S2 is relevant, while when A < (blue-detuned laser), the 
second condition is always satisfied and only S\ matters. 



B. Correlation matrix of the quantum fluctuations 
of the system 

The steady state of the bipartite quantum system 
formed by the vibrational mode of interest and the fluc- 
tuations of the intracavity mode can be fully charac- 
terized. In fact, the quantum noises £ and a m are 
zero-mean quantum Gaussian noises and the dynamics 
is linearized, and as a consequence, the steady state 
of the system is a zero-mean bipartite Gaussian state, 
fully characterized by its 4 x 4 correlation matrix (CM) 
Vij = ((ui(oc)uj(oa) + Uj(oo)ui(oa))) /2. Starting from 
Eq. ijHJ), this steady state CM can be determined in two 
equivalent ways. Using the Fourier transforms Uj(cj) of 
Ui(t), one has 



Vij(t) 



u(t) = Au(t) + n(t), 



(8) 



dtodu}' 

47T 



Ui(u))Uj(u) ) + Uj(u>)Ui(u))) 

(12) 
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Then, by Fourier transforming Eq. (JSj) an d the correlation 
functions of the noises, Eqs. (fj| and ([5]), one gets 



(Ui(uj)Uj(uj') + Uj{<jj')Ui(uj)) 



M{uj)D{u)M(J) t 5(<J + u'), 

ij 



where we have defined the 4x4 matrices 

v-l 



M(u) = (ilu + Ay 



and 



D(u) 











°\ 





^coth(^) 














K 





Vo 








J 



(13) 
(14) 

(15) 
(16) 



The S(oj + oj') factor is a consequence of the stationarity 
of the noises, which implies the stationarity of the CM 
V: in fact, inserting Eq. (fl"3]) into Eq. (fl2|) . one gets that 
V is time-independent and can be written as 



V = <IujM{uj)D{uj)M{uj) ] . 



(17) 



It is however reasonable to simplify this exact expression 
for the steady state CM, by appropriately approximating 
the thermal noise contribution -D22 i n Eq. (|16p . In 
fact hsT/h ~ 10 11 s _1 even at cryogenic temperatures 
and it is therefore much larger than all the other typical 
frequency scales, which are at most of the order of 10 9 Hz. 
The integrand in Eq. (|17[) goes rapidly to zero at u> ~ 10 11 
Hz, and therefore one can safely neglect the frequency 
dependence of Dxi{^>) by approximating it with its zero- 
frequency value 



y m u ( Huj 
coth 1 



\2k B T 



2k B T 



7m (2n + l), (18) 



where n = (exp{huj m / k B T} — 1) is the mean thermal 
excitation number of the resonator. 

It is easy to verify that assuming a frequency- 
independent diffusion matrix D is equivalent to make 
the following Markovian approximation on the quantum 
Brownian noise 

WW) + £(*')£(*)> /2 * j m (2n + l)S(t - f ), (19) 

which is known to be valid also in the limit of a very high 
mechanical quality factor Q = u) m /y m —> 00 (68|. Within 
this Markovian approximation, the above frequency do- 
main treatment is equivalent to the time domain deriva- 
tion considered in [53] which, starting from the formal 
solution of Eq. ([8]) , arrives at 

/•OO POO 

V^(co)=V/ da I ds'M^M^Duis-s 1 ), 
k,i Jo Jo 

(20) 



where D kl (s — s') = ((n k (s)ni(s') + rn(s')n k (s))) /2 is 
the matrix of the stationary noise correlation functions. 
The Markovian approximation of the thermal noise on 
the mechanical resonator yields Dki{s — s r ) = Dki5(s — 
s'), with D = Diag[0,7 m (2n+ 1), k,k], so that Eq. ([20]) 
becomes 



V = I dsM(s)DM(s) T , 
Iq 



(21) 



which is equivalent to Eq. (JT7J) whenever D does not de- 
pend upon oj. When the stability conditions are satisfied 
(M(oo) = 0), Eq. (f2"Tj) is equivalent to the following Lya- 
punov equation for the steady-state CM, 



AV + VA T = -D, 



(22) 



which is a linear equation for V and can be straightfor- 
wardly solved, but the general exact expression is too 
cumbersome and will not be reported here. 



III. OPTOMECHANICAL ENTANGLEMENT 
WITH THE INTRACAVITY MODE 

In order to establish the conditions under which the 
optical mode and the mirror vibrational mode are entan- 
gled we consider the logarithmic negativity EV, which 
can be defined as 169] 



where rj 
with S(V) 



£jv = max[0, — In 2rj ] , 



2 -i/2 



(23) 



1 1/2' 



1/2 



E(V) - [S(V) 2 -4dety]" 

det V m + det V c — 2 det V mc , and we have 
used the 2x2 block form of the CM 



V 



V 1 



mc 

v c 



(24) 



Therefore, a Gaussian state is entangled if and only if 
?/~ < 1/2, which is equivalent to Simon's necessary 
and sufficient entanglement non-positive partial trans- 
pose criterion for Gaussian states [70j , which can be writ- 
ten as 4dety < S - 1/4. 



A. Correspondence with the down-conversion 
process 

As already shown in [28|, [2^, [3(| many features of the 
radiation pressure interaction in the cavity can be under- 
stood by considering that the driving laser light is scat- 
tered by the vibrating cavity boundary mostly at the first 
Stokes (ujq — uj m ) and anti-Stokes (u)q + u) m ) sidebands. 
Therefore we expect that the optomechanical interaction 
and eventually entanglement will be enhanced when the 
cavity is resonant with one of the two sidebands, i.e., 
when A = ±w„, . 
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It is useful to introduce the mechanical annihilation 
operator 5b = (Sq + idp)/y/2, obeying the following QLE 



.G 



5b = -iuj m 5b -J^.(5b- Stf) + (5a f + 5a) + -±= 

(25) 

Moving to another interaction picture by introducing the 
slowly-moving tilded operators Sb(t) — Sb(t)e~ Wmt and 
da(t) = 5a(t)e~ lAt , we obtain from the linearized version 
of Eq. (3b) and Eq. (25]) the following QLEs 

Sb = -2^L (sb-slU 2 ^ 

2 V 

^gt^A-h^t + SSLel (u m -A)t^ (26) 



—n5a 



G 



(27) 



Note that we have introduced two noise operators: i) 
a m {t) = a m (i)e lAt , possessing the same correlation func- 
tion as a m (t); ii) b m (t) = £{t)e iUmt / '^2 which, in the 
limit of large uj m , acquires the correlation functions (7lj 



(b m '\t)b m (t')) = nS(t-t'), 
(b m {t)b m ^(t')) = [n + l]5(t-t'). 



(28) 
(29) 



27)) are still equivalent to the linearized QLEs 
but now we particularize them by choosing 
If the cavity is resonant with the Stokes side- 



Eqs. (26] 
of Eq. § 

A = ±uj, 

band of the driving laser, A 

5l = _2™«56 + I^^t e 2^ mi 



m , one gets 
■ G ,~t 



5a = -nSa + i—Stf + i^6be 2luJmt 



2nd 1 



(30) 
(31) 



while when the cavity is resonant with the anti-Stokes 
sideband of the driving laser, A = u) n 



5b 



5 a 



7rc 



■5b 



■5~tfe 2tuJmt 



2 

—n5a + - 



i — 5b^e 
2 



j m , one gets 
■x-Sa 



2na l 



(32) 
,(33) 



From Eqs. (3"0")) - (f3"l"j) we see that, for a blue-detuned driv- 



ing laser, A = 



the cavity mode and mechanical 



resonator are coupled via two kinds of interactions: i) a 
down-conversion process characterized by Sb^da^ + 5a5b, 
which is resonant and ii) a beam-splitter-like process 
characterized by 5w5a + Sa^Sb, which is off resonant. 
Since the beam splitter interaction is not able to entan- 
gle modes starting from classical input states (72j, and 
it is also off-resonant in this case, one can invoke the 
rotating wave approximation (RWA) (which is justified 



in the limit of oj m 3> G, k) and simplify the interaction 
to a down conversion process, which is known to gen- 
erate bipartite entanglement. In the red-detuned driv- 
ing laser case, Eqs. (J3"2"]) - (|3"3")) show that the two modes 
are strongly coupled by a beam-splitter-like interaction, 
while the down-conversion process is off-resonant. If one 
chose to make the RWA in this case, one would be left 
with an effective beam splitter interaction which cannot 
entangle. Therefore, in the RWA limit ui m 3> G, k, the 
best regime for strong optomechanical entanglement is 
when the laser is blue-detuned from the cavity resonance 
and down-conversion is enhanced. However, as it will be 
seen in the following section, this is hindered by instabil- 
ity and one is rather forced to work in the opposite regime 
of a red-detuned laser where instability takes place only 
at large values of G. 



B. Entanglement in the blue-detuned regime 

The CM of the Gaussian steady state of the bipar- 
tite system, can be obtained from Eqs. (f5U|) - (f5T]) and 
Eqs. Eqs. (3S |) -(33 )) in the RWA limit, with the techniques 
of the former section (see also [73| ) 



V = V ± 









± 

14 





















(34) 



where the upper (lower) sign corresponds to the blue- 
(red-)detuned case, and 



V, 



:s:i 



1 2G 2 k[1/2± (n+ 1/2)] 

2 + (7m + 2k) {2 lm n T G 2 ) ' 
G 2 lm [n + 1/2 ± 1/2] 
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2 ( 7m + 2k) {2 lm K T G 2 ) ' 
2Gj m n[n + l/2±l/2] 
(j m + 2 K ) (2j m K =p G 2 ) ' 



(35a) 
(35b) 
(35c) 



For clarity we have included the red-detuned case in the 
RWA approximation and we see that det V^ c = T(V^) 2 , 
i.e., is non- negative in this latter case, which is a suf- 
ficient condition for the separability of bipartite states 
[70l ]. Of course, this is expected, since it is just the beam- 
splitter interaction that generates this CM. Thus, in the 
weak optomechanical coupling regime of the RWA limit, 
entanglement is obtained only for a blue-detuned laser, 
A = — uj m . However, the amount of achievable optome- 
chanical entanglement at the steady state is seriously lim- 
ited by the stability condition of Eq. (TTTa ). which in the 
RWA limit A = —u> m k, 7m , simplifies to G < y/2nj m . 
Since one needs small mechanical dissipation rate 7 m in 
order to see quantum effects, this means a very low max- 
imum value for G. The logarithmic negativity E^f is an 
increasing function of the effective optomechanical cou- 
pling G (as expected) and therefore the stability condi- 
tion puts a strong upper bound also on Ej^. In fact, it is 
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possible to prove that the following bound on Em exists 
Em < In 



1 + G/V2^ 



1 + n 



(36) 



showing that Em < In 2 and above all that entanglement 
is extremely fragile with respect to temperature in the 
RWA limit because, due to the stability condition, Em 
vanishes as soon as n > 1. 



C. Entanglement in the red-detuned regime 

We conclude that, due to instability, one can find sig- 
nificant optomechanical entanglement, which is also ro- 
bust against temperature, only far from the RWA regime, 
in the strong coupling regime in the region with posi- 
tive A, because Eq. (jllb ) allows for higher values of G. 
This is confirmed by Fig. [2j where Em is plotted versus 
the normalized detuning A/ui m and the normalized input 
power P/Po, (Po — 50 mW) at a fixed value of the cavity 
finesse F = Fq = 1.67 x 10 4 in (a), and versus the normal- 
ized finesse F/Fq and normalized input power P/Pq at a 
fixed cavity detuning A = uj m in (b). We have assumed 
an experimentally achievable situation, i.e., a mechanical 
mode with oj m /2ir = 10 MHz, Q = 10 5 , mass m = 10 ng, 
and a cavity of length L = 1 mm, driven by a laser with 
wavelength 810 nm, yielding Go = 0.95 kHz and a cavity 
bandwidth k = 0.9uj m when F = Fq. We have assumed 
a reservoir temperature for the mirror T = 0.4 K, corre- 
sponding to h ~ 833. Fig. [2^ shows that Em is peaked 
around A ~ u m , even though the peak shifts to larger 
values of A at larger input powers P. For increasing P at 
fixed A, Em increases, even though at the same time the 
instability region (where the plot suddenly interrupts) 
widens. In Fig.[2p we have fixed the detuning at A = uo m 
(i.e., the cavity is resonant with the anti-Stokes sideband 
of the laser) and varied both the input power and the 
cavity finesse. We see again that Em is maximum just 
at the instability threshold and also that, once that the 
finesse has reached a sufficiently large value, F ~ Fq, Em 
roughly saturates at larger values of F. That is, one gets 
an optimal optomechanical entanglement when n ~ uo m 
and moving into the well-resolved sideband limit k <C Lu m 
does not improve the value of Em ■ The parameter region 
analyzed is analogous to that considered in [53| , where it 
has been shown that this optomechanical entanglement 
is extremely robust with respect to the temperature of 
the reservoir of the mirror, since it persists to more than 
20 K. 



D. Relationship between entanglement and cooling 

As discussed in detail in [H, [2iJ [3(| HH the same 
cavity-mechanical resonator system can be used for re- 
alizing cavity-mediated optical cooling of the mechanical 
resonator via the back-action of the cavity mode [23| . In 
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FIG. 2: (Color online) (a) Logarithmic negativity Em versus 
the normalized detuning A/ui m and normalized input power 
P/Po, (Po = 50 mW) at a fixed value of the cavity finesse 
F = F Q = 1.67 x 10 4 ; (b) Em versus the normalized finesse 
F/Fo and normalized input power P/Pq at a fixed detuning 
A = uj m . Parameter values are ui m /2-K = 10 MHz, Q — 10 5 , 
mass m — 10 ng, a cavity of length L = 1 mm driven by 
a laser with wavelength 810 nm, yielding Go = 0.95 KHz 
and a cavity bandwidth k = 0.9u> m when F — Fq. We have 
assumed a reservoir temperature for the mirror T = 0.4 K, 
corresponding to n ~ 833. The sudden drop to zero of Em 
corresponds to entering the instability region. 



particular, back-action cooling is optimized just in the 
same regime where A ~ u> m . This fact is easily explained 
by taking into account the scattering of the laser light 
by the oscillating mirror into the Stokes and anti-Stokes 
sidebands. The generation of an anti-Stokes photon takes 
away a vibrational phonon and is responsible for cooling, 
while the generation of a Stokes photon heats the mirror 
by producing an extra phonon. If the cavity is resonant 
with the anti-Stokes sideband, cooling prevails and one 
has a positive net laser cooling rate given by the differ- 
ence of the scattering rates. 

It is therefore interesting to discuss the relation be- 
tween optimal optomechanical entanglement and opti- 
mal cooling of the mechanical resonator. This can eas- 
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ily performed because the steady state CM V deter- 
mines also the resonator energy, since the effective sta- 
tionary excitation number of the resonator is given by 
n eff = (V n + V22 - 1) /2 (see Ref. [H for the exact ex- 
pression of these matrix elements giving the steady state 
position and momentum resonator variances). In Fig. [3] 
we have plotted n e ff under exactly the same parameter 
conditions of Fig. O We see that ground state cooling is 
approached (n e ff < 1 ) simultaneously with a significant 
entanglement. This shows that a significant back-action 
cooling of the resonator by the cavity mode is an impor- 
tant condition for achieving an entangled steady state 
which is robust against the effects of the resonator ther- 
mal bath. 

Nonetheless, entanglement and cooling are different 
phenomena and optimizing one does not generally op- 
timize also the other. This can be seen by comparing 
Figs. [2] and [3j E_\f is maximized always just at the in- 
stability threshold, i.e., for the maximum possible op- 
tomechanical coupling, while this is not true for n e //, 
which is instead minimized quite far from the instability 
threshold. For a more clear understanding we make use 
of some of the results obtained for ground state cooling 
in Refs. [H, [29|, [3(|. In the perturbative limit where 
G <C uj m ,K, one can define scattering rates into the 
Stokes (A+) and anti-Stokes (A^) sidebands as 



A± = 



G 2 k/2 



so that the net laser cooling rate is given by 



T = A_ - A+ > 0. 



(37) 



(38) 




The final occupancy of the mirror mode is consequently 
given by [H,[li,[33 



FIG. 3: (Color online) (a) Effective stationary excitation 
number of the resonator n e ff versus the normalized detun- 
ing A/com and normalized input power P/Po, (Po = 50 mW) 
at a fixed value of the cavity finesse F = Fo — 1.67 x 10 4 ; 
(b) n e ff versus the normalized finesse F/Fo and normalized 
input power P/Pq at a fixed detuning A = ui m . Parameter 
values are the same as in Fig. [2] Again, the sudden drop to 
zero corresponds to entering the instability region. 



A+ 



7m + r 7 m + r 



(39) 



IV. OPTOMECHANICAL ENTANGLEMENT 
WITH CAVITY OUTPUT MODES 



where the first term in the right hand side of the above 
equation is the minimized thermal noise, that can be 
made vanishingly small provided that j m -C T, while 
the second term shows residual heating produced by 
Stokes scattering off the vibrational ground state. When 
T 7m^j the lower bound for n e ff is practically set by 
the ratio A + /T. However, as soon as G is increased for 
improving the entanglement generation, scattering into 
higher order sidebands takes place, with rates propor- 
tional to higher powers of G. As a consequence, even 
though the effective thermal noise is still close to zero, 
residual scattering off the ground state takes place at a 
rate that can be much higher than A + . This can be seen 
more clearly in the exact expression of (Sq 2 ) — Vn given 
in [30], which is shown to diverge at the threshold given 
byEq. CETb). 



The above analysis of the entanglement between the 
mechanical mode of interest and the intracavity mode 
provides a detailed description of the internal dynamics of 
the system, but it is not of direct use for practical appli- 
cations. In fact, one typically does not have direct access 
to the intracavity field, but one detects and manipulates 
only the cavity output field. For example, for any quan- 
tum communication application, it is much more impor- 
tant to analyze the entanglement of the mechanical mode 
with the optical cavity output, i.e., how the intracavity en- 
tanglement is transferred to the output field. Moreover, 
considering the output field provides further options. In 
fact, by means of spectral filters, one can always select 
many different traveling output modes originating from 
a single intracavity mode and this gives the opportunity 
to easily produce and manipulate a multipartite system, 
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eventually possessing multipartite entanglement. 



A. General definition of cavity output modes 

The intracavity field Sa(t) and its output are related 
by the usual input-output relation [63| 



t (t) = V2^6a(t)-a in (t), 



(40) 



where the output field possesses the same correlation 
functions of the optical input field a ln (t) and the same 
commutation relation, i.e., the only nonzero commutator 
is [a out (t),a out (t')i] =6{t-t'). From the continuous out- 
put field a out (t) one can extract many independent opti- 
cal modes, by selecting different time intervals or equiv- 
alently, different frequency intervals (see e.g. [7J])- One 
can define a generic set of N output modes by means of 
the corresponding annihilation operators 



dsg k (t- S )a out (s), k = l,...N, (41) 



where g k (s) is the causal filter function defining the fc-th 
output mode. These annihilation operators describe N 
independent optical modes when \a° ut {t), a° k ut (t)'~\ = 6j k , 
which is verified when 



dsgj(s)*g k {s) = S jk , 



(42) 



i.e., the N filter functions g k (t) form an orthonormal set 
of square-integrable functions in [0, oo). The situation 
can be equivalently described in the frequency domain: 
taking the Fourier transform of Eq. (|4Tj) . one has 



dt 



a° k ut {t)e iut 



2K~g k {u J )a out {uj), 



where g k {uj) is the Fourier transform of the filter func- 
tion. An explicit example of an orthonormal set of filter 
functions is given by 



0(t)-6(t-r) _ iQkt 
9k(t) = 7= e 



(44) 



(8 denotes the Heavyside step function) provided that £l k 
and r satisfy the condition 



r^ 7 - — fit 



2tt 



-p, integer p. 



(45) 



These functions describe a set of independent optical 
modes, each centered around the frequency Q k and with 
time duration r, i.e., frequency bandwidth ~ 1/t, since 



9k(v) 



(w - n k )r/2 



(46) 



When the central frequencies differ by an integer multiple 
of 27t/t, the corresponding modes are independent due 
to the destructive interference of the oscillating parts of 
the spectrum. 

B. Stationary correlation matrix of output modes 

The entanglement between the output modes defined 
above and the mechanical mode is fully determined by 
the corresponding (2N + 2) x (2iV + 2) CM, which is 
defined by 



where 



'(*) = \ {uT\t)uf\t) + u^{t)u1 ut {t)) , (47) 



(48) 



= (5q(i),5p(t), Xr'(i), Yr\t), . . . , X°^{t),YZ vt {t)) T 

is the vector formed by the mechanical position 
and momentum fluctuations and by the ampli- 
tude (X° ut (t) = [a° k ut (t)+a° k ut (ty] /V2), and phase 



'(i)*] /iV2) quadratures of the 



N output modes. The vector u out (t) properly describes 
N + 1 independent CV bosonic modes, and in particu- 
lar the mechanical resonator is independent of (i.e., it 
commutes with) the N optical output modes because 
the latter depend upon the output field at former times 
only (s < t). From the definition of u out (t), of the out- 
put modes of Eq. (|41[) , and the input-output relation of 
Eq. ([40]) one can write 



dsT lk (t-s)u e k xt (s) 
dsT lk (t- S )nt xt ( S ), (49) 



where 



u ext {t) = (6q(t), 6p(t),X(t), Y(t), X(t),Y(t)f 

(50) 

is the 2N + 2-dimcnsional vector obtained by extending 
the four-dimensional vector u(t) of the preceding section 
by repeating N times the components related to the op- 
tical cavity mode, and 



1 



'2k 



(0,0, X m (t),Y m (t), . . .,X in (t), Y m (t)) J 



(51) 

is the analogous extension of the noise vector n{t) of the 
former section without however the noise acting on the 
mechanical mode. In Eq. (|4"9"]) we have also introduced 
the (2N + 2) x (27V + 2) block-matrix consisting of N + 1 
two-dimensional blocks 
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T(t) = 



(8{t) 
5(t) 







2nRegi{t) 
2K,Imgi(t) 







_0 
2Klm(7i (t) 
2K,Regi(t) 










2KReg 2 {t) 
2nlvag2{t) 








2Klmg2 (t) 
2KReg 2 {t) 



(52) 



r 



Using Fourier transforms, and the correlation function of 
the noises, one can derive the following general expression 



J 



-y-out 



dwT(uj) 



M ext (uj) + 



2k 



for the stationary output correlation matrix, which is the 
counterpart of the 4x4 intracavity relation of Eq. (fl~7|) 



D ext {uj) 



P„ 



2k 



(53) 



where P out — Diag[0, 0, 1, 1, . . .] is the projector onto 
the 2iV-dimensional space associated with the output 
quadratures, and we have introduced the extensions cor- 
responding to the matrices M(u>) and D(u>) of the former 
section, 



with 



and 





M 


ext {oj) 


= (iuj 


+ A exty 


-l 




(54) 


/ 













. 


"\ 




— w m 




G 


G 


. 












— K 


A 


. 








G 





-A 


-AC 


. 




■ (55) 













-ft 


A . 






G 








-A 


— K . 






V 












■■) 





D ext (uj) 







^f coth (^ 










V 







K K 

K K 

K K 

K K 



(56) 

A deeper understanding of the general expression for 
yout £ (|53p is obtained by multiplying the terms 
in the integral: one gets 



V ouf = duf(Lo)M ext (co)D ext (Lo)M ext (Lo^f{Lu)^ + ^ + - <kof{u) M ext {Lo)R out +R out M ext {uj)^ T{lo)\ (57) 



r 



where R out = P out D ext (uj) / k = D ext (to) P out / k and we 
have used the fact that 



^f(co)P out D ext (co)P out f(u;^ 
Ak* 



(58) 



the second noise term and eventually, even increase the 
optomechanical entanglement with respect to the intra- 
cavity case. We shall analyze this fact in the following 
section. 



The first integral term in Eq. (|57p is the contribution 
coming from the interaction between the mechanical res- 
onator and the intracavity field. The second term gives 
the noise added by the optical input noise to each output 
mode. The third term gives the contribution of the cor- 
relations between the intracavity mode and the optical 
input field, which may cancel the destructive effects of 



C. A single output mode 

Let us first consider the case when we select and de- 
tect only one mode at the cavity output. Just to fix 
the ideas, we choose the mode specified by the filter 
function of Eqs. (|4"4"|) and (I46p . with central frequency 
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fl and bandwidth r _1 . Straightforward choices for this 
output mode are a mode centered either at the cavity 
frequency, fl — ui c — loq, or at the driving laser frequency, 
f2 = (we are in the rotating frame and therefore all 
frequencies are referred to the laser frequency ujq), and 
with a bandwidth of the order of the cavity bandwidth 
t _1 ~ k. However, as discussed above, the motion of the 
mechanical resonator generates Stokes and anti-Stokes 
motional sidebands, consequently modifying the cavity 
output spectrum. Therefore it may be nontrivial to de- 
termine which is the optimal frequency bandwidth of the 
output field which carries most of the optomechanical en- 
tanglement generated within the cavity. The cavity out- 
put spectrum associated with the photon number fluctua- 
tions S(oj) = (Sa(ujy 8a(u>)) is shown in Fig.[4j where we 
have considered a parameter regime close to that con- 
sidered for the intracavity case, i.e., an oscillator with 
ujm/2-K = 10 MHz, Q = 10 5 , mass m — 50 ng, a cavity 
of length L = 1 mm with finesse F — 2 x 10 4 , detuning 
A = u> m , driven by a laser with input power P = 30 
mW and wavelength 810 nm, yielding Go = 0.43 kHz, 
G = 0.41u> m , and a cavity bandwidth k = 0.75w m . We 
have again assumed a reservoir temperature for the mir- 
ror T — 0.4 K, corresponding to n ~ 833. This regime 
is not far but does not corresponds to the best intra- 
cavity optomechanical entanglement regime discussed in 
Sec. IIIII In fact, optomechanical entanglement monoton- 
ically increases with the coupling G and is maximum just 
at the bistability threshold, which however is not a conve- 
nient operating point. We have chosen instead a smaller 
input power and a larger mass, implying a smaller value 
of G and an operating point not too close to threshold. 
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FIG. 4: (Color online) Cavity output spectrum in the case of 
an oscillator with uj m /2ir — 10 MHz, Q = 10 5 , mass m = 50 
ng, a cavity of length L = 1 mm with finesse F = 2 x 10 4 , 
detuning A = uj m , driven by a laser with input power P = 
30 mW and wavelength 810 nm, yielding Go = 0.43 kHz, 
G = 0.41oj m , and a cavity bandwidth « = 0.75w m . We have 
again assumed a reservoir temperature for the mirror T = 0.4 
K, corresponding to n ~ 833. In this regime photons are 
scattered only at the two first motional sidebands, at loq ±oj m . 

In order to determine the output optical mode which 
is better entangled with the mechanical resonator, we 



study the logarithmic negativity Ej*j associated with the 
output CM V out of Eq. (JF2) (for N = 1) as a function 
of the central frequency of the mode £1 and its band- 
width r _1 , at the considered operating point. The re- 
sults are shown in Fig. [5j where EV is plotted versus 
Q/u> m at five different values of e = ruj m . If e < 1, i.e., 
the bandwidth of the detected mode is larger than u> m , 
the detector does not resolve the motional sidebands, and 
Em has a value (roughly equal to that of the intracavity 
case) which does not essentially depend upon the cen- 
tral frequency. For smaller bandwidths (larger e), the 
sidebands are resolved by the detection and the role of 
the central frequency becomes important. In particular 
Ejj becomes highly peaked around the Stokes sideband 
Q = —Lu m , showing that the optomechanical entangle- 
ment generated within the cavity is mostly carried by 
this lower frequency sideband. What is relevant is that 
the optomechanical entanglement of the output mode is 
significantly larger than its intracavity counterpart and 
achieves its maximum value at the optimal value e ~ 10, 
i.e., a detection bandwidth t" 1 ~ u> m /10. This means 
that in practice, by appropriately filtering the output 
light, one realizes an effective entanglement distillation 
because the selected output mode is more entangled with 
the mechanical resonator than the intracavity field. 




FIG. 5: (Color online) Logarithmic negativity Ej^ of the CV 
bipartite system formed by mechanical mode and a single cav- 
ity output mode versus the central frequency of the detected 
output mode Q./io m at five different values of its inverse band- 
width e — u> m T. The other parameters are the same as in 
Fig. [4] When the bandwidth is not too large, the mechani- 
cal mode is significantly entangled only with the first Stokes 
sideband at luq — u) m . 

The fact that the output mode which is most entangled 
with the mechanical resonator is the one centered around 
the Stokes sideband is also consistent with the physics of 
two previous models analyzed in Refs. [55|,[75|]. In [75[ an 
atomic ensemble is inserted within the Fabry-Perot cavity 
studied here, and one gets a system showing robust tri- 
partite (atom-mirror-cavity) entanglement at the steady 
state only when the atoms are resonant with the Stokes 
sideband of the laser. In particular, the atomic ensemble 
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and the mechanical resonator become entangled under 
this resonance condition, and this is possible only if en- 
tanglement is carried by the Stokes sideband because the 
two parties are only indirectly coupled through the cav- 
ity mode. In [55| . a free-space optomechanical model is 
discussed, where the entanglement between a vibrational 
mode of a perfectly reflecting micro-mirror and the two 
first motional sidebands of an intense laser beam shined 
on the mirror is analyzed. Also in that case, the mechan- 
ical mode is entangled only with the Stokes mode and it 
is not entangled with the anti-Stokes sideband. 

By looking at the output spectrum of Fig. [3J one can 
also understand why the output mode optimally entan- 
gled with the mechanical mode has a finite bandwidth 
t _1 ~ w m /10 (for the chosen operating point). In fact, 
the optimal situation is achieved when the detected out- 
put mode overlaps as best as possible with the Stokes 
peak in the spectrum, and therefore r _1 coincides with 
the width of the Stokes peak. This width is determined 
by the effective damping rate of the mechanical resonator, 
7^/^ = 7 m + r, given by the sum of the intrinsic damping 
rate j m and the net laser cooling rate T of Eq. (|38l It is 
possible to check that, with the chosen parameter values, 
the condition e = 10 corresponds to r _1 ~ 7^. 

It is finally important to analyze the robustness of 
the present optomechanical entanglement with respect to 
temperature. As discussed above and shown in [53j |. the 
entanglement of the resonator with the intracavity mode 
is very robust. It is important to see if this robustness 
is kept also by the optomechanical entanglement of the 
output mode. This is shown by Fig. [51 where the entan- 
glement £V °f the output mode centered at the Stokes 
sideband £1 = — oj m is plotted versus the temperature of 
the reservoir at two different values of the bandwidth, 
the optimal one e = 10, and at a larger bandwidth e = 2. 
We see the expected decay of Ej^ for increasing tempera- 
ture, but above all that also this output optomechanical 
entanglement is robust against temperature because it 
persists even above liquid He temperatures, at least in 
the case of the optimal detection bandwidth e = 10. 



D. Two output modes 

Let us now consider the case where we detect at the 
output two independent, well resolved, optical output 
modes. We use again the step-like filter functions of 
Eqs. (I44[) and (f^o]) , assuming the same bandwidth t -1 
for both modes and two different central frequencies, f2i 
and SI2 , satisfying the orthogonality condition of Eq. (Hoi) 
Ql — £1% = 2p7TT~ 1 for some integer p, in order to have 
two independent optical modes. It is interesting to an- 
alyze the stationary state of the resulting tripartite CV 
system formed by the two output modes and the mechan- 
ical mode, in order to see if and when it is able to show 
i) purely optical bipartite entanglement between the two 
output modes; ii) fully tripartite optomechanical entan- 
glement. 
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FIG. 6: (Color online) Logarithmic negativity Ej^ of the CV 
bipartite system formed by mechanical mode and the cavity 
output mode centered around the Stokes sideband f2 = — uo m 
versus temperature for two different values of its inverse band- 
width e = u> m T. The other parameters are the same as in 
Fig.H 



The generation of two entangled light beams by means 
of the radiation pressure interaction of these fields with a 
mechanical element has been already considered in var- 
ious configurations. In Ref. [76| . and more recently in 
Ref. [12] , two modes of a Fabry-Perot cavity system with 
a movable mirror, each driven by an intense laser, are 
entangled at the output due to their common pondero- 
motive interaction with the movable mirror (the scheme 
has been then generalized to many driven modes in [77| ) . 
In the single mirror free-space model of Ref. [55|, the 
two first motional sidebands are also robustly entangled 
by the radiation pressure interaction as in a two-mode 
squeezed state produced by a non-degenerate parametric 
amplifier (78[. Robust two- mode squeezing of a bimodal 
cavity system can be similarly produced if the movable 
mirror is replaced by a single ion trapped within the cav- 
ity [ZH. 

The situation considered here is significantly differ- 
ent from that of Refs. [H, M lli, 

which require 

many driven cavity modes, each associated with the cor- 
responding output mode. In the present case instead, 
the different output modes originate from the same sin- 
gle driven cavity mode, and therefore it is much simpler 
from an experimental point of view. The present scheme 
can be considered as a sort of "cavity version" of the 
free-space case of Ref. [55[ , where the reflecting mirror is 
driven by a single intense laser. Therefore, as in [55lF7c|. 
one expects to find a parameter region where the two out- 
put modes centered around the two motional sidebands of 
the laser are entangled. This expectation is clearly con- 
firmed by Fig. where the logarithmic negativity Ej^f 
associated with the bipartite system formed by the out- 
put mode centered at the Stokes sideband (fii = —u> m ) 
and a second output mode with the same inverse band- 
width (e = Lo m T — IO71") and a variable central frequency 
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f2, is plotted versus Q/uj m . E^ is calculated from the 
CM of Eq. d57J) (for N = 2), eliminating the first two 
rows associated with the mechanical mode, and assum- 
ing the same parameters considered in the former sub- 
section for the single output mode case. One can clearly 
see that bipartite entanglement between the two cavity 
outputs exists only in a narrow frequency interval around 
the anti-Stokes sideband, Q = u> mi where £V achieves its 
maximum. This shows that, as in [H|,[78[, the two cav- 
ity output modes corresponding to the Stokes and anti- 
Stokes sidebands of the driving laser are significantly en- 
tangled by their common interaction with the mechanical 
resonator. The advantage of the present cavity scheme 
with respect to the free-space case of [EH [lH is that the 
parameter regime for reaching radiation-pressure medi- 
ated optical entanglement is much more promising from 
an experimental point of view because it requires less 
input power and a not too large mechanical quality fac- 
tor of the resonator. In Fig. [51 the dependence of E^ 
of the two output modes centered at the two sidebands 
fi = ±w m upon their inverse bandwidth s is studied. We 
see that, differently from optomechanical entanglement 
of the former subsection, the logarithmic negativity of 
the two sidebands always increases for decreasing band- 
width, and it achieves a significant value (~ 1), compa- 
rable to that achievable with parametric oscillators, for 
very narrow bandwidths. This fact can be understood 
from the fact that quantum correlations between the two 
sidebands are established by the coherent scattering of 
the cavity photons by the oscillator, and that the quan- 
tum coherence between the two scattering processes is 
maximal for output photons with frequencies luq i oj m . 
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FIG. 8: (Color online) Logarithmic negativity Eu of the bi- 
partite system formed by the two output modes centered at 
the Stokes and anti-Stokes sidebands {Q. = ±oj m ) versus the 
inverse bandwidth e = w m r. The other parameters are the 
same as in Fig. [4] 



by plotting, for the same parameter regime of Fig. [8j 
Ejy versus the temperature T at two different values of 
the inverse bandwidth (e = 107r, 1007r). We see that 
this purely optical CV entanglement is extremely robust 
against temperature, especially in the limit of small de- 
tection bandwidth, showing that the effective coupling 
provided by radiation pressure can be strong enough 
to render optomechanical devices with high-quality res- 
onator a possible alternative to parametric oscillators for 
the generation of entangled light beams for CV quantum 
communication. 
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FIG. 7: (Color online) Logarithmic negativity Ex of the bi- 
partite system formed by the output mode centered at the 
Stokes sideband (Sli = — uj m ) and a second output mode with 
the same inverse bandwidth (e = ui m T — 10tv) and a vari- 
able central frequency Q, plotted versus (l/u m . The other 
parameters are the same as in Fig. [4] Optical entanglement 
is present only when the second output mode overlaps with 
the anti-Stokes sideband. 

In Fig. [5] we analyze the robustness of the entangle- 
ment between the Stokes and anti-Stokes sidebands with 
respect to the temperature of the mechanical resonator, 



1. 



0.5 



1 1 










1 1 










T 1 

£ = \0x 

1_ | 










£ = 

1 1 


100* 









\ 





















, 






20 



T(K) 



40 



60 



FIG. 9: (Color online) Logarithmic negativity Ej^ of the two 
output modes centered at the Stokes and anti-Stokes side- 
bands (fi = ±w m ) versus the temperature of the resonator 
reservoir, at two different values of the inverse bandwidth 
e = u> m T. The other parameters are the same as in Fig. [4] 

Since in Figs. [7] and [5] we have used the same parameter 
values for the cavity-resonator system used in Fig. [5l we 
have that in this parameter regime, the output mode cen- 
tered around the Stokes sideband mode shows bipartite 
entanglement simultaneously with the mechanical mode 



13 



and with the anti-Stokes sideband mode. This fact sug- 
gests that, in this parameter region, the CV tripartite 
system formed by the output Stokes and anti-Stokes side- 
bands and the mechanical resonator mode could be char- 
acterized by a fully tripartite-entangled stationary state. 
This is confirmed by Fig. I10( where we have applied the 
classification criterion of Ref . [8(| , providing a necessary 
and sufficient criterion for the determination of the en- 
tanglement class in the case of tripartite CV Gaussian 
states, which is directly computable in terms of the eigen- 
values of appropriate test matrices These eigenval- 
ues are plotted in Fig. [TU] versus the inverse bandwidth 
e at A = uj m in the left plot, and versus the cavity de- 
tuning A/tu m at the fixed inverse bandwidth e = tt in 
the right plot (the other parameters are again those of 
Fig. 2|) . We see that all the eigenvalues are negative in 
a wide interval of detunings and detection bandwidth of 
the output modes, showing, as expected, that we have a 
fully tripartite-entangled steady state. 



A 





FIG. 10: (Color online) Analysis of tripartite entanglement. 
The minimum eigenvalues after partial transposition with re- 
spect to the Stokes mode (blue line), anti-Stokes mode (green 
line) and the mechanical mode (red line) are plotted versus 
the inverse bandwidth e at A = u) m in the left plot, and ver- 
sus the cavity detuning A /ui m at the fixed inverse bandwidth 
e — 7r in the right plot. The other parameters are the same 
as in Fig. 2] These eigenvalues are all negative in the studied 
intervals, showing that one has fully tripartite-entanglement. 



Therefore, if we consider the system formed by the 
two cavity output fields centered around the two mo- 
tional sidebands at ujq ± uj m and the mechanical res- 
onator, we find that the entanglement properties of its 
steady state are identical to those of the analogous tri- 
partite optomechanical free-space system of Ref. [55| . In 
fact, the Stokes output mode shows bipartite entangle- 
ment both with the mechanical mode and with the anti- 
Stokes mode, the anti-Stokes mode is not entangled with 
the mechanical mode, but the whole system is in a fully 
tripartite-entangled state for a wide parameter regime. 
What is important is that in the present cavity scheme, 
such a parameter regime is much easier to achieve with 
respect to that of the free-space case. 



V. CONCLUSIONS 



We have studied in detail the entanglement proper- 
ties of the steady state of a driven optical cavity coupled 
by radiation pressure to a micromechanical oscillator, 
extending in various directions the results of Ref. [53| . 
We have first analyzed the intracavity steady state and 
shown that the cavity mode and the mechanical clement 
can be entangled in a robust way against temperature. 
We have also investigated the relationship between entan- 
glement and cooling of the resonator by the back-action 
of the cavity mode, which has been already demonstrated 
recently in Refs. [ill ElJlMlJlMljM and discussed 
theoretically in Refs. ^lllliniSllallir We have seen 
that a significant back-action cooling is a sufficient but 
not necessary condition for achieving entanglement. In 
fact, intracavity entanglement is possible also in the op- 
posite regime of negative detunings A where the cav- 
ity mode drives and does not cool the resonator, even 
though it is not robust against temperature in this lat- 
ter case. Moreover, entanglement is not optimal when 
cooling is optimal, because the logarithmic negativity is 
maximized close to the stability threshold of the system, 
where instead cooling is not achieved. 



We have then extended our analysis to the cavity out- 
put, which is more important from a practical point of 
view because any quantum communication application 
involves the manipulation of traveling optical fields. We 
have developed a general theory showing how it is pos- 
sible to define and evaluate the entanglement properties 
of the multipartite system formed by the mechanical res- 
onator and N independent output modes of the cavity 
field. 

We have then applied this theory and have seen that in 
the parameter regime corresponding to a significant in- 
tracavity entanglement, the tripartite system formed by 
the mechanical element and the two output modes cen- 
tered at the first Stokes and anti-Stokes sideband of the 
driving laser (where the cavity output noise spectrum is 
concentrated) shows robust fully tripartite entanglement. 
In particular, the Stokes output mode is strongly entan- 
gled with the mechanical mode and shows a sort of en- 
tanglement distillation because its logarithmic negativity 
is significantly larger than the intracavity one when its 
bandwidth is appropriately chosen. 

In the same parameter regime, the Stokes and anti- 
Stokes sideband modes are robustly entangled, and the 
achievable entanglement in the limit of a very narrow de- 
tection bandwidth is comparable to that generated by a 
parametric oscillators. These results make the present 
cavity optomechanical system very promising for the re- 
alization of CV quantum information interfaces and net- 
works. 
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